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Abstract 

In this paper we study a model for phase segregation consisting in a sistem of a partial 
and an ordinary differential equation. By a careful definition of maximal solution to the 
latter equation, this system reduces to an Allen-Cahn equation with a memory term. 
Global existence and uniqueness of a smooth solution are proven and a characterization of 
the u-hmit set is given. 
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1 Problem setting 

The Allen-Cahn equation 

Kdtp-Ap + f'ip) = (1.1) 

is meant to describe evolutionary processes in a two-phase material body, including phase seg- 
regation: p, with p{x, t) e [0, 1], is an order-parameter field interpreted as the scaled volumetric 
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density of one of the two phases, k > is a mobility coefficient, and / denotes a double-well 
potential confined in (0, 1) and singular at cndpoints. The derivation of this equation offered 
by Gurtin in [3] is based on a balance of contact and distance microforces: 

div^ + 7r + 7 = (f.2) 

and on a 'purely mechanical' dissipation inequality restricting the free-energy growth: 

dti^<w, w:=-TTdtp + ^-W{dtp), (1.3) 

where the distance niicroforcc is split in an internal part tt and an external part 7, ^ denotes the 
microscopic stress vector, and w specifies the (distance and contact) internal microworkingQ 
the Coleman-Noll compatibility of the constitutive choices 

7r = ^(p,Vp,atp), ^-i(p,Vp,9tp), and ^ = V^(p, Vp) = /(p) + ^jVpp (1.4) 

with the dissipation inequality (jl.3p yields 

n{p, Vp, dtp) = -f'{p) - k(p, Vp, dtp)dtp, ^(p, Vp, dtp) = Vp, (1.5) 

and hence the Allen-Cahn equation (jl.ip . for k(p, Vp, dtp) — k and 7 = 0. 

One of us proposed in [7] a modified version of Gurtin's derivation, where the dissipation 
inequality p.Sp is dropped and the microforce balance ()1.2p is coupled with the microenergy 
balance 

dte — e + w, e :— — div h + a, (1.6) 
and the microentropy imbalance 

dtri > ~ div h + a, h ph, a := pa. (1-7) 

The salient new feature of this approach to phase-segregation modeling is that the microentropy 
inflow {h,(7) is deemed proportional to the microenergy inflow {h,a) through the chemical 
potential p, a positive field; consistently, the free energy is defined to be 

i;:=e-p-\, (1.8) 

with the chemical potential playing the same role as coldness in the deduction of the heat 
equation^ Combination of p.6p - p.8p gives: 

9tV < -vip-^J+ p-'h ■Vp-ndtp + ^- y{dtp), (1.9) 

an inequality that replaces for (|1.3p in restricting constitutive choices that can now be more 
general than those in (|1.4p . On taking all of the constitutive mappings delivering tt, ^, rj, and h 
depending on the list p, Vp, dtp, and the chemical potential p, and on choosing 

ij ^ ^{p,Vp,p) = -pp+f{p) + i|Vpp, (1.10) 

^In [2] the balance of microforces is stated under form of a principle of virtual power for microscopic motions. 

^Just as absolute temperature can be seen as a macroscopic measure of microscopic agitation, its inverse - 
the coldness - measures microscopic quiet; likewise, the chemical potential can be seen as a macroscopic measure 
of microscopic organization. 
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compatibility with (jl.9|) yields 

tt{p, Vp, dtp, m) = p- f'{p) - k{p, Vp, dtp)dtp, ^{p, Vp, ^tp, ^) = Vp, 
'7(P, Vp, 9tp, ^) = -p p, /i(p, Vp, 9tp, ^) = 0. 

With the use of p. lip and, once again, of the additional constitutive assumptions that the 
mobility k is a positive constant and the external distance microforcc 7 is null, the microforce 
balance ()1.2p and the energy balance p.6p become, respectively, 

Kdtp-Ap + f'ip)^p (1.12) 

and 

dti-p.^p)^p{nidtpy+a). (1.13) 

This nonlinear system consists of a parabolic PDE and a first-ordcr-in-timc ODE and is to 
be solved for the order-parameter field p and the chemical potential field p; formally, setting 
p = restitutes the standard AUen-Cahn equation (jl.ip . We supplement system (|1.12p - (|1.13p 
with the homogeneous Neumann condition 

dnp = on the body's boundary (1-14) 

(here dn denotes the outward normal derivative) and with the initial conditions 

p\t=Q = Po bounded away from , p|t=o = Po > . (1-15) 

Note that, in view of the third of relations (jl.lip . the microcntropy cannot exceed the level 
from below anywhere at any time, and that the corresponding prescribed initial field 

v\t=o ^ Vo ^ -f^Po (1-16) 

is nonpositive- valued. 

Remark 1.1. The last of p. lip implies that both the energy influx and the entropy influx are 
everywhere null for all times. This result, that has a pivotal role in reducing the energy balance 
to an ODE, is a direct consequence of assuming that the energy-influx mapping h, just as all 
the other constitutive mappings, be independent of the gradient of the chemical potential. This 
assumption, that precludes energy and entropy diffusion, does not seem appropriate in the case 
of higher-order phase segregation models of Cahn-Hilliard type, like the one proposed in [7]. 

2 Solution strategy and summary of contents 

The aim of our paper is a mathematical investigation of problem (|1.12p - p.l5p . The key idea 
of our strategy is to attack the problem sequentially, the ODE flrst, then the PDE. 

To do so, we introduce a change of variable that is expedient to give (|1.13p plus (|1.16p the 
form of a parametric initial- value problem. The change of variable in question is: 



(2.1) 
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whence 
it leads to 



(2.2) 



a.e + ^^^^v/e = 0, ^Uo = Co, (2.3) 

a Cauchy problem for ^(x, •) parameterized on the space variable x and on the field p(x, •). 
The general form of this problem is discussed in the Appendix. Suffice it to note here that 
(|2.3p exhibits the Peano phenomenon and has infinitely many solutions; among them, we pick 
a suitably defined maximal solution ^ (or y/£^, see definition (|3.6p - (|3.7p ). having the desirable 
property to stay positive as long as is possible. Next, we transform (|1.12p into 

Kdtp-Ap + f'ip)-^^^0, (2.4) 

that is, an AUen-Cahn equation for p(x, ■) with the additional term —\f^Jp\ since the factor 
is implicitly defined in terms of p as the maximal solution of (j2.3p , ()2.4p must be regarded as 
an integrodifferential equation. We prove existence, regularity and uniqueness of the solution to 
(j2.4p subject to the boundary condition (|1.14p and the initial condition (jl.lSP i by means of a 
fixed-point argument, taking advantage of the iterated Contraction Mapping Principle. Crucial 
to success is that dtp he a, priori uniformly bounded in the space-time domain; we show that 
this is the case by applying standard regularity arguments for parabolic equations. 

For a detailed discussion of the problem transformation sketched here above, and for a precise 
exposition of our main results, we refer the reader to Section [31 Our well-posedness results are 
proved in Section [H Section E] is devoted to an investigation of the long-time behavior of the 
solution; we prove that uniquely converges to some function ip^o and that any element p^o 
of the limit set solves the stationary problem 

Kdtp- Apoo + f'{Poo)-Voo^^ = 0, (2.5) 



supplemented by suitable homogeneous Neumann boundary conditions. 



3 Main results 

We begin by specifying once and for all the class of data we consider; further assumptions of 
local importance will be stated when needed. Our problem is formulated over a space-time 
cylinder 

QT^nx[0,T) with T e (0,-t-oo), (3.1) 

where il is an open, bounded and connected set ofM.^ {N > 1). with a smooth boundary P (we 
use the notation Qt :~ ^ x [0,t) for every t G (0, +oo]). As to the coarse-grain free energy /, 
we split it as follows: 

< / = /i + /2, where /i, /s : (0, 1) ^ R are C^-functions, (3.2) 
/i is convex, is bounded, lim f'{r) = — oo, and lim f'{r) = +oo. (3.3) 

r\0 ryi 
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We regard /2 as a smooth perturbation of the singular convex part /i of /, which is well 
exemplified by 

/i(r) =rlnr+(l-r)ln(l-r) for r G (0, 1) . (3.4) 
As to the energy source and the initial data, we assume that 

aeL^iQr), po, Co e L°°(f7), < po < 1 and > a.e. in ^2. (3.5) 

Finally, we recall that the mobility k is a given positive constant. 

With this, we take up the forward Cauchy problem (|2.3p . Clearly. ^ must be nonnegative. 
We notice that, if we look for a strictly positive ^ (for given p > and > 0), the Cauchy 
problem ()2.3p has a unique local solution. On the contrary, uniqueness is no longer guaranteed 
if we allow ^ to be just nonnegative. On the other hand, every nonnegative local solution can 
be extended to a global solution. Therefore, we select a (global) solution to problem (|2.3|1 
according to the following maximality criterion (for a justification, see the Appendix): 



\/^{x, t) = sup {w{x, t) : ui G §*((T, p)} for (x, <) G Qt, where (3.6) 
S*(CT,eo,p) {we W^'''(0,r;L\l7)) : u;(0) = v/^, a.e. in Qt, 

dtw -^{nidtpf +a)/{2p^''^) a.e. where w > o|. (3.7) 

Accordingly, the maximal ^ satisfies: 



V'?(a^,i) = / a*{x,s)ds, (3.8) 

Jo 

where 

f i^\dtp{x, s)p +a(x,s) 

*/ N , it t(a;,s)>0, ,„ 

a*{x,s) := <^ 2^^{^ ' ' (3.9) 

[ otherwise. 

At this point, we replace fu, by \f^Jp in p.l2p and supplement the equation (|2.4p we get with 
the boundary and initial conditions for p given by, respectively, (|1.14p and the first of (|1.15|) . Of 
the so-obtained initial/boundary value problem we give the following variational formulation: 

for 

V:=H^{^) and E:^l?{^), (3.10) 

seek a field p such that: 

pG ifi(0,r;iJ)nC"([0,T];y); (3.11) 

p(0) = po, 0<p<l a.e. in Qt, - + ^ ^°°(Qt) ; (3.12) 

p 1 - p 

K [ dtp{t)z+ f Vp(t)-Vz+ / r{p{t))z- [ (e(t)/p(t))'/'z-0 (3.13) 
Jn Jn Jn Jn 

for a.a. t G (0,T), for every z G V, and for ^ given by 

Remark 3.1. We regard the initial/boundary value problem p.lip - p.l3p as an essentially 
intcgrodifferential AUen-Cahn equation in the sole unknown p. We note, in particular, that 
()3.13p has a precise meaning, because S,^/"^ G L^{Qt) and p^^^'^ G L°°{Qt) (at least) whenever 
p satisfies (|3.1ip and a G L'^{Qt)- 
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Here is our first result (the symbol ( • ) denotes the negative part). 
Theorem 3.2. Assume that: 

Sei°°(Qoo) and (J- e L\0,oo;L°^\n)); — + G L°°(f7), (3.14) 

Po 1 - Po 

po e H^{n), dnPQ = onT, and Apo e L°°{n). (3.15) 
Then, for every T G (0, +oo), problem (|3.1ip - (|3.13p has a unique solution. Furthermore, 

pe LP{0,T;W^-'P{n)) for every p < +O0, dtpeL^iQr), and ^eL°°{QT). (3.16) 
Finally, there exist constants p*,p* G (0, 1) and ^* > such that 

Pt: < p < p* and £, < S,* a.e. in Qt; (3-17) 
these constants can be chosen independently of T . 

Remark 3.3. Let A/q be an upper bound for Assume that := inf is strictly positive, 
and take M such that |i9tp| < M a.e. in Qt (sec p.l6p ). Then, relation (|6.10p in the Appendix 
implies that ^(i) > at least for t < 2y/^^/{HiM'^ + A/q). In such a case, the pair (p, ^) yields 
a (local) solution (p, p, rf) to the original problem in a strong sense. 

Our second result concerns the long-time behavior of the solution p to problem (|3.1ip - (|3.13p : 
it ensures that the elements of the w-limit of every trajectory are steady states. To state this 
result properly, we have to describe the stationary problem associated to p.lip - (|3.13p . We 
introduce ipoo : ^ [0, +00), by means of the following formula: 

Voo{x) lim v?OM0 foi' 2: G il, where ^ is given by p.6p - p.9p (3.18) 

(our next theorem ensures that such a limit actually exists, under form of a bounded function 
on ri). The stationary problem consists in finding 

Poo G y and < Poo < P* a-e. in Q. (3.19) 

such that 

/ Vpoo • Vz+ / /'(poc)z- / ^z = for every zGV. (3.20) 

Theorem 3.4. Under the assumptions of Theorem \3.2l let p be the unique global solution to 
problem p.lip - p.l3p . Then, the limit ()3.18p exists for a. a. a; G andipoo G £°°(r2). Moreover, 
the Lo -limit defined by 

uj{p) '■= {p°° G H : p°° = lim p{tn) strongly in H for some {tn} /* +00} (3.21) 

n — 'OO 

is non-empty, compact, and connected in the strong topology of H . Finally, every element 
p°° G uj{p) coincides with a solution poo to the stationary problem (j3.19p - (|3.20p . 

Remark 3.5. One can wonder whether /i can be a more general potential from Convex 
Analysis. Actually this is the case, since we may replace the monotone part /{ of /' by a 
graph a. Precisely, we may assume that 

a is a maximal monotone graph in M x M, (3.22) 
with £'(0;) = (0, 1), limr\o ck°('') = "OO, and limryi Q!*^(r) = +cxd, (3.23) 



7 



where D(a) stands for the domam of a and a^[r) denotes the element of a{r) having mmimuni 
modulus for r e (0,1) (see, e.g., [I] p. 28]). Accordingly, /i is replaced by a convex l.s.c. 
function S : M ^ (— oo, +oo] such that da — a, so that f a + f2 > 0- 

When thinking of such a generalization, we have to introduce a selection C of a{p) with 
some regularity and we have to come up with a convenient replacement for the variational 
equation (|3.13p . Here is a suitably general formulation: 

find C&L^iQr) and C G "(p) 

a.c. in Qt, such that (3.24) 
Jn Jn Jn Jn 

for a.a. t e (0, T), for every v eV, and for ^ given by ([SJ]) . (3.25) 

An analogous modification is due for the stationary problem (|3.19p - (|3.20p . that may be replaced 
by the following problem: 

find Coo G L^(il) and Coo S a(poo) a.e. in il, such that (3.26) 

/ Vpoo • V2 + / (Coo + /2(Poo)) z - f ^ z = for every z g V. (3.27) 
Jn Jn Jn vPoo 

With such measures. Theorem 13.21 can be extended as far as existence and regularity are con- 
cerned. Precisely, we can prove that there is a global solution (p, C)i with C G L°"{Q oo) and p 
satisfying the same regularity requirements and bounds as in Theorem 13.21 However, wc can- 
not prove uniqueness. Furthermore, Theorem 13.41 holds for every (possibly nonunique) global 
solution to the generalized problem satisfying the same bounds as above. We sketch how to 
achieve such generalizations in the forthcoming Remarks 14.91 and 15.31 



4 Existence and uniqueness 

In this section, wc prove Theorem 13. 2 1 This is a rather complicated task, that we now delineate. 
First of all, wc show that every solution satisfies the last part of the statement, i.e., that some 
kind of maximum principle holds. Then, wc show that we can count a priori on more regularity 
than that specified in (j3.16p . Finally, by looking for solutions satisfying such stronger properties, 
only, wc prove existence and uniqueness using a fixed point argument. In the preliminary steps, 
we find convenient to deal with an auxiliary problem. In the whole section, it is understood 
that the assumptions of Theorem 13.21 are satisfied. 

Construction of the crucial constants. We find the constants p*, p* , and noting that 
our procedure actually yields values that do not depend on T, as stated in the last part of 
Theorem 13.21 For convenience, we set 

M2:= sup |/^(r)|, (4.1) 

re(0,l) 

and choose G (0, 1) and Co — have 

/{(p*)<-^^2, Po>P* and Co<Co a.e. inf7, (4.2) 

due to (|3.3p . p.5p . and p.l4p . Moreover, on accounting for the second of (|3.14p . we define 
C* > as follows: 

\f¥ ■■= \/^+7ri=\\^^\\LHo,oc;L'^(n))- (4.3) 
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Finally, using the last of (j3.3p and (j3.5|) once more, we choose p* E (0, 1) such that 

fUp*)- Ah and po<p* a.c. in fi. (4.4) 



At this point, we are ready to prove the last part of Theorem l3.2l At the same time, with a 
view towards the fixed point argument we are going to use later on, we prepare some auxiliary 
material. So, we show that 

p* < p < p* and ^ < ^* a.e. in Qt (4.5) 

for any solution p to the variational equation (|3.13p . but with (|3.6P " P.9p replaced by something 
else. 

The auxiliary problem. From the previous section, we see that y/£^, rather than ^, plays 
the main role. Hence, wc define $ : !?($) W'^'^{0, T; L^^)) ^ L'^iQr) as follows. We set: 

!)($):= {v e i?^(0,T;i/), v>0 a.e. in Qt, f/w G L°°(Qt)} (4.6) 

and, for v € £*(*&), we denote by <^{v) the function given by 

tp{x,t) = s\xp{w{x,t) : w e 3'iv)} for {x,t) € Qt, (4.7) 

where we have set 

3^{v) := eW^'\0,T;L\n)) : w{0) = ^/^, w>0 a.e. in Qt, 

ajw==-(K(9tw)2 + CT)/(2i;i/2) a.e. where u; > o|. (4.8) 

By arguing as in the Appendix, we see that ip := $(w) is the square root of the maximal solution 
to the Cauchy problem: 

dt^^-[{^idtv)' + a)/i2v'/^))^ and ^(0) = ?o, 

and it is characterized by 

ifiix^t) = y/^o{x) - / a*{x,s)ds, (4.9) 
Jo 

where 

n*(j^^) '^\9tv{x, s)\^ + ct {x _^ p{x,t)>0, a*(x,s):=0 otherwise. (4.10) 

2^/v{x, s) 

Note that ip actually belongs to W^-^{0,T; L^{fl)) n L°°{Qt)- Then, the auxiliary problem is 
obtained as follows. For a given v e £*($), we require that p satisfies: 

pei/i(0,r;iJ)nC"([0,T];V^), p(0) = po, (4.11) 
< p< 1 a.e. in Qt, f'{u) G L^(Qt), and p^^/^ g L^{Qrj.)- (4.12) 

K f dtp{t) z + / Vp{t) ■Vz+ f f'{p{t)) z - f z (4.13) 



for a. a. t E (0,T) and every z gV, where (p — $(u). 

Therefore, problem (j3.1ip - (j3.13p is equivalent to the auxiliary problem, provided v ~ p and 
p ~ Y^, and provided that some stronger regularity requirements arc granted. 
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Lemma 4.1. Let v G D{^) and assume that p satisfies (j4.11|) - (|4.14p . Then, we have that 

p> p^, a.e. in Qt- (4-14) 
In particular, this is true if v is any solution p to problem (|3.1ip - (j3.13p . 

Proof. The proof wc give is quite standard. Let g : M ^ R be Lipschitz continuous, nonde- 
creasing, and such that g{r) < for r < p:, and g{r) ~ for r > p»; furthermore, let G be the 
primitive of g that vanishes at p*. Now, we write ()4.14|) at t = s and test it by z := g{p{s)). 
Then, we integrate over (0,t) with respect to s, where t e (0,T) is arbitrary. We have: 

^ / G{pit)) + [ Vp- Vg{p) + / if'ip) - gip) [ G{p^). (4.15) 

Jn JQt JQt JO 



The right-hand side vanishes by (|4.2[) and the first two terms on the left-hand side are non- 
negative. As to the third, we can replace Qt by its subset where p < p^,. Due to (|3.2p - (|3.3p 
and (|4.ip - (|4.2p . a.e. in such a subset we have that: 

/'(p)-^</'(p)</{(p*) + M2<0 and g{p)<0, (4.16) 

whence the nonnegativity of the third integral in (|4.15p . We conclude that G{p{t)) = a.e. 
in n for every t € [0, T], i.e., that p > p* a.e. in Qt- □ 

Lemma 4.2. Assume v G Di^) and v > p*. Then, we have that 

^{v) < a.e. m Qt- (4.17) 

In particular, ^ < ^* a.e. in Qt for every solution (p, ^) to problem (j3.1ip - (|3.13p . 

Proof. We set := $(t;) for brevity and stipulate that 

X is the characteristic function of the subset of Qt where > 0. (4-18) 
Then, dUD-gini) yield 



2v{s) 







2 Jo \/v[S) ^y/P* 

In particular, ^ < ^* if w = p and (p, solves problem p.lip - p.l3p . because the inequality 
p > p^t has been already proved for every solution. □ 

Lemma 4.3. Assume v S D{^) and v > p* a.e. in Qt and let p satisfy ()4.1ip - ()4.14p . Then, 
we have that 

p < p* a.e. in Qt. (4-19) 
In particular, this is true if v is any solution p to problem p. lip - (|3.13p . 
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Proof. Let g : M ^ K. be Lipschitz continuous, nondecreasing, such that g{r) > for r > p* 
and g{r) = for r < p*; and let G be the primitive of g that vanishes at p*. Then, (|4.15p holds 
with the new g and G; once again, the only trouble conies from the third term on the left-hand 
side. However, (a.e.) in the subset of Qt where p > p*, wc have that 



g{p)>0 and f'{p)-^>f[(p)-^ + f^{p)>f[(p*)-^-M,_>0, (4.20) 

thanks to the estimate (|4.17p and the definition of p* given by (|4.4p . We conclude that G{p{t)) = 
a.e. in il for every t G [0,T], i.e., that p < p* a.e. in Qt- □ 



Remark 4.4. We note that the last sentence of Theorem 13.21 is completely proved. 

Lemma 4.5. Let v £ £'($) be such that p* < v < p* a.e. in Qt, and let p satisfy ()4.1ip - ()4.14p . 

Then, we have that 

\\dtp\\LP{QT) ^ fo"^ ewer?/ p e (1, +oo), (4.21) 
IIpIIlp(o,T;W2.p(o)) < R'p for every pe (l,+oo), (4.22) 

where the constants Rp and R'^ depend only on the structure of our problem, the initial data pQ 
and ^Q, and p. In particular, all this is true if v is any solution p to problem p.lip - p.l3p . 

Proof. Observe that 

/ dtp{t)z + / Vp(t) • Vz = / F{t)z and p{0) = po 
Jn Jn Jn 

for a.a. t £ (0, T), for every zeV, and for F Lp/^ - f'{p). Owing to gTl]), gH]), gTH), 
and to our assumptions on / and poi we have that 

|F|<^+ sup |/'(r)| a.e. in Qt and ||po||L=^(n) < P*- (4.23) 

Therefore, we can apply the general L^'-rcgularity theory (see, e.g., Thm. 9.1, p. 341]) and 
deduce that (|4.2ip - (|4.22[l hold. Moreover, thanks to the above lemmas, the first of ()3.16p is 
proved. □ 

Lemma 4.6. Let v e -D($) be such that p* < v < p* a.e. in Qt, and let p satisfy (|4.11[) - (|4.14p . 
Moreover, assume that 

\\dtv\\Lp(QT) < Rp for every p G (1, +oo), (4.24) 
with the same Rp as in ()4.2ip . Then, 

\\dtp\\L-iQ,)<Roo, (4.25) 

for some constant Rrx, depending only on the structure of our problem and the initial data po 
and ^Q. In particular, this is true if v is any solution p to problem (j3.1ip ~ (|3.13p . 
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Proof. Here we use the stronger regularity assumption (j3.15|) . We proceed formally, because the 
argument that would make the calculation rigorous is quite standard. We differentiate (|4.14p 
with respect to time and obtain: 

K / dtu{t)z+ / Vu(i)-Vz+ / u{t)z= i F{t)z for a.a. t G (0, T) and every z G (4.26) 
Jo, Jn Jn Jn 

where we have set set u := dtp and 

F := dtp - f"{p)dtp - \^p'^'^dtp -\x {K\dtv? + a) {vp)-''\ 

with X as in ()4.18p . Now, the initial value of u is known explicitely through the differential 
equation for p, and is given by the formula: 

Ku{Q) = Apo - f'iPo) + VIJVPo ■ 

Hence, w(0) belongs to L°°{Q) and its L°°-norm is estimated by a known constant, due to p.lSp . 
Now, we observe that the following estimate holds (see (SJ Thm. 7.1, p. 181]): 

Ml^(q^) < Cgmax{||u(0)|U»(o) + (4.27) 

where q and Cq are as follows: q is any real number satisfying q > rj^/{rN ~ 1), where rpf = 
max{2, 1 + (A^/2)}; Cq is a constant depending only on and q. Hence, we are led to find an L''- 
estimate of F for some q satisfying the above inequality. Actually, we can choose anyone of such 
values of q and get an estimate using the previous lemma and assumptions p.l4p and ()4.24p . 
If we denote by C a suitable constant that could be computed in terms of our assumptions on 
the problem structure and the constants p*, p*, and we have indeed that 

WFWlHQt) < C iWdtphHQT) + \\idtv)^\\L'>{QT) + W^hoiQT)) 
< C{Rq + i?2g + I|0'!Il9(Qt))' 

for any q E (1, +oo). Hence, (|4.27|) provides a value of Roo satisfying (|4.25p . □ 
Remark 4.7. We note, in particular, that now all the conditions listed in p.l6p are proved. 

Towards the fixed point argument. Due to the above lemmas, we can confine ourselves to 
look for solutions satisfying all the bounds we have proved, i.e., belonging to the set X defined 
here below; 

X:^{veV: p.<v<p* a.c. in Qt, \\dtv\\LPiQ^) < Rp ^ v <^ (1, +oo]}, (4.28) 

where the constants involved are of course the same as before, and where 

V := ff\0,T;i7)nC"([0,r];F). (4.29) 

We regard X as a metric subspace of the Banach space V, i.e., we consider the metric d on X 
defined by 

(f{u,v):=[ \dtu ^ dtv\^ + sup \\u{t) - v{t)\\lr foru,z;eX. (4.30) 
Jqt te[o,T] 
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Thus, X is complete: each one of the conditions specified in the definition (|4.28p of X defines 
a closed subset of V, and such is X, being the intersection of a family of closed sets. Next, for 
f S X, we observe that v G -0(3"), and we define the mapping : X ^ X as follows: 

\E'(u) is the unique solution p to problem (|4.1ip - (|4.14p . (4-31) 

It is clear that, for every u g X, the auxiliary problem (|4.1ip - (|4.14p has a unique solution p. 
Indeed, if we set Lp := $(«), we sec that, for a. a. (x^t) € Qt, the nonlinear function r > 
/'(r) — if(x,t)/y/r, defined for p* < r < p*, is Lipschitz continuous; moreover, ip is bounded. 
Hence, the standard theory for regular parabolic equations yields existence and uniqueness and 
we conclude that '^{v) is well defined. Furthermore, the precise choice of all constants in (|4.28p 
exactly ensures that '^{v) belongs to X, due to the previous lemmas, whence ^ actually maps 
X into itself. Finally, it is clear that an element of X is a solution to (j3.1ip - (j3.13p if and only 
if it is a fixed point for ^ . Hence, it suffices to prove that '^^ o ■ ■ ■ o [k times) is a 

contraction for k large enough. The rest of the proof is devoted to do that. A key point of our 
argument is the following lemma. 

Lemma 4.8. For i ~ 1,2, pick ai G Lj^^.^, +oo), and let yi be the maximal solution to the 
Cauchy problem (|6.ip with a = Oi. Then, we have 



sup |\/2/i(s) - Vy^is)] < / \ai{s) - a2{s)\ dr for every t > 0. (4.32) 
se[o,t] JO 



Proof. We first prove that 

^ - Vy^l < |ai - 02! a.e. in (0, +00). (4.33) 

To this end, we choose everywhere defined representatives of a,;, set Lpi := ^Jyl for brevity, 
and notice that the functions tpi and \(pi — (p2\ are locally absolutely continuous on [0, +00). 
Therefore, there exists an exceptional set E C [0, +00) having zero Lebesgue measure, such 
that the derivatives at t of the above functions exist for every t S (0, +00) \ E. Moreover, 
we have that </5-(i) = —ai{t) if ipi{t) > 0. Let us fix a point t outside of E and prove (|4.33p 
at t. We distinguish three cases. In the first one, we have (pi{t) > for i = 1,2. Then, 
\^i ~ <^2| = i(¥'i — (^2) in a neighborhood of t and ipl{t) = ~ai{t) for i = 1, 2, whence 

Wi - 'P2\'{t) = ±{^i - ip2Y{t) = T(ai - a2){t) < \ai{t) - a2(i)|. 

Assume now (pi{t) = for i = 1,2. Then, (p[{t) — (see (|6.7p ') and the desired inequality trivially 
follows. In the last case, we have, e.g., ipi{t) > = (p2{t). We derive that (p[{t) = —ai{t). 
On the other hand, as 'fi2{t) = 0, we see that (|6.8p implies 02 (t) > 0. Therefore, noting that 
(pi — ^2 > in SL neighborhood of t, we deduce that 

Ifi - V2\'{t) = ((^1 - ^2Y{t) = ^\{t) - v'2{t) < -ai(t) + a2{t) < |ai(t) - a2{t)\, 

and ()4.33p is completely proved. Now, we derive (|4.32p . We fix t > and s £ [0,t]. Then, 
(1433]) yields 

Wviis) - Vy^is)] = / \Vyi - Vy^\' (r) dr < / |ai(r) - a2(r)| dr < / |ai(T) - a2(T)| dr, 
Jo Jo Jo 



and (|4.32p immediately follows. 



□ 
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Conclusion. We take vi,V2 G X and set for convenience (pi := ^{vi) and pi := '^{vi) for 
i = 1,2. Tlien, we write the equality in (|4.14p for v = vi and test the difference by dt{pi — P2)- 
Then, we integrate over (0, t) for an arbitrary t e (0, T) and add the same integral to both sides 
for convenience. If we set p := pi — p2 and use a similar notation for Lp and w, we obtain: 

^1 \dtP? + \\\p{t)\\l=l {p-f'{pi) + riP2))dtp+f Up^'/^^p2P2'^'')dtpiAM) 
JQt ^ JQt JQt ^ ' 

We now pass to estimate the right-hand side of the last relation. In order to simplify the 
notation, we use the same symbol c for different constants (even in the same formula) that 
only depend on the problem structure, the data, and T; denotes those such constants that 
depend, in addition, on the parameter i5 e (0, 1). As the functions / and ( • are Lipschitz 

continuous on the interval [p.^, p*], and as estimate (|4.17|) holds for both we have that 



(p-/'(pi)+/'(P2))aiP + [p^^p-^l^ ^p2P2^l'')dtP 

<cf \p\\dtp\+ f pi\p\\dtp\+ f P2'^'\p\\dtp\ 

JQt JQt JQt 

<S [ \dtp\'+cs I \p\' + cs I Iv^P, 



for every 5 G (0, 1). By combining with (j4.34p and choosing 5 small enough, we easily deduce 
that ^ ^ 

\dtP? + \\p{t)fv<c I \\p{s)\\]jds + c( Ms)\\lds. 

Jo JO 

On the other hand, we can compare p.6p - p.9p . (|6.5p - (|6.7p . and (|4.6p - (|4.10p and apply Lem- 
ma HiHl Thus, for every s g (0,t) and a.e. in fi, we have that 



\tpis)\ = \(pi{s) - ip2is)\ 

dT 



1 

< - 

1 

< - 
~ 2 



+ f i ^2'^'i^)\\dMrr~\dtV2{ 



<c {l + \a{T)\)\v{T)\dT + c \dtv{T)\dT<c \v{T)\dT + c \dtv{T)\ dr, 
Jo Jo JO JO 

since < Vi < p* . Hence, 

WfisWu < c||w||i2(Q^) + c\\dtv\\l2,Q^ <c \dtif + c\\v\\cono,s];V)- 

JQs 

Therefore, we deduce that 

^ \dtp? + \\pml<cj\p{s)\\lds + cj*^(^j^ |5tif + ds, 
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and the Gronwall lemma easily yields: 

\dtp{' + l|pilco([o.«];y) 1^*^'!' + ll'^llco([o,.];y)) for every t € [0, T], 

with a precise constant C*. This means that 

(A)! / (/ l^*^""' " "'^'^ ^ """^ " ^2|lcO([0,.];V)) '^S: 



< 



for every t £ [0,T] with k = 1. Arguing by induction on k, it is straightforward to prove that 
the above inequality holds for every fc > 1. We conclude that 

(f{-^''{vi) ~ ^^(v2)) < — (f(vi,V2) for every integer A: > 1, (4.35) 

k\ 

whence, 'I'*^ is a contraction on X for k large enough. This completes the proof. 

Remark 4.9. Here, we breafly sketch how to modify the above proof in order to achieve the 
generalization mentioned in Remark 13.51 where ()3.13|) was replaced by (|3.24p - (j3.25p . Accord- 
ingly, we consider the auxiliary problem obtained by assuming (|3.24p for C and modifying (j4.14p 
as follows: 

K / dtp{t)z+ f Vp(t)-Vz+ / (C + /2(pW))^- / (eW/pW)'^'^ = (4.36) 
for a. a. t £ (0, T) and every z £ V, where ^ is given by (|4.7p . 



However, it is convenient to consider approximating problems as well. Precisely, we write a^i^p) 
instead of ( in equations (|3.25p and (|4.37p . where a,, denotes the Yosida regularization of a 
at level e € (0, 1) (see, e.g. [U p. 28]). We will speak of the e-gencralized problem and of the 
e-auxiliary problem, respectively. Now, we briefly show how to obtain first uniform bounds as 
in the previous lemmas, and then existence, for the generalized problem. As far as the choice 
of the crucial constants is concerned, in conditions (|4.2p and (|4.4p we replace /{ and M2 by q;° 
and M2 + 1, respectively. 

As |ae(?')| < |Q!°(r)|, and as ae(r) converges to a°(r) when e tends to zero for every r € (0, 1), 
we see that the inequalities 

aeip*)<-M2 and a,{p*) - ^ > M2 



hold true for e small enough. Therefore, it is easy to see that Lemmas I4.1H4.3I still hold for 
each e-problem, i.e., that the a priori bounds (j4.5p are fulfilled in the new situation. 



As to Lemma [4.51 we observe that 

a" (r' )< a£(r) < «"(/') whenever < r' < < r < p* < r" < 1. (4.37) 

Therefore, bounds (|4.5I) imply uniform bounds for C, = ae{p) that can play the role of ()4.23p . 
Thus, we see that (|4.2ip holds for the solution to each of the e-problems uniformly with re- 
spect to e. 
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As far as the generalization of Lemma 14.61 is concerned, time differentiation is allowed for 
e-problcms. However, we replace (|4.26p by the following equality 

K / dtu{t)z+ / Vu(i)-Vz+ / u{t)z+ / a'^{p)u{t)z^ / F{t)z, 
Jq Jq Jn Jn Jn 

and we modify the previous definition of F by writing in place of /". In other words, we 
replace the term — /"(/o)(?tP on the right-hand side by the term a'^{p)u on the left-hand side. 
As a possible technique for (|4.27p consists in testing the above equation by z := Z{u), where 
Z is monotone and vanishes at in order to get recursive L^'-estimates by a Moser's argument, 
and as the integral involving a'^ gives a nonnegative contribution in such a procedure, the type- 
()4.25|) bound one finds for the solutions to our e-problems is uniform with respect to e. All 
this ensures that the definition of X can be done in the present case, and it is independent 
of e. However, the constant C* we find in applying the fixed point argument docs depend 
on s. Nevertheless, this is enough to conclude for the existence of a unique solution to the 
approximating problem among the functions p belonging to X. Since uniqueness among all 
solutions can be proved for every fixed e > 0, we see that our argument constructs a global 
solution defined in the whole of [0, 4-oo). Moreover, for every fixed T, p^ satisfies a number 
of a priori estimates uniformly with respect to e, and (|3.17p and (j4.37p imply that a uniform 
X°°-estimate holds for := a^ipe)- Therefore, modulo standard arguments, we see that (pe, Ce) 
sequentially converges in the proper topology to a pair (p, () and that (p, C) is a global solution 
of the generalized problem. Unfortunately, in the new situation the previous uniqueness proof 
does not work, because ()3.25p cannot be differentiated with respect to time. 



5 Long-time behavior 

In this section, we prove Theorem 13.41 Hence, we choose data satisfying the prescribed condi- 
tions and pick the corresponding unique global solution (p, ^). Our proof is organized as follows. 
The first lemma establishes the first assertion of the theorem, i.e., that the function ipo^ given 
by (|3.18p is well-defined and bounded. In the next one, we derive some new a priori estimates, 
which ensure that the cj-limit uj^p) given by (j3.2ip has the desired properties. Finally, we 
conclude by the announced characterization of iu{p). We find it convenient to set: 

(fi := and X := the characteristic function of the subset of Qoo where (fi > 0, (5.1) 

and we notice that the function tpoo to be studied is the pointwise limit of ip as time goes to 
infinity. Furthermore, we recall that estimates (|4.14p and (|4.19p hold for p. 

Lemma 5.1. The limit p.lSp is well defined and Lp^o is bounded. 
Proof. Equations (|i^ - (|i?TUl) and definitions ((F?T|) yield: 

Jo 2p(s) 
where we have set ^ 

A_(x,t) f X(a;, s) ^ ^^1^ ds 

Jo 2^p(x, s) 



16 



and 

K\dtp{x,s)\'^ + a+{x,s) 



\+{x,t) := X{x,s) ■ — ' ■ — ds 

Jo 2y^ p{x,s) 

for a. a. {x,t) £ Qoc- To prove the assertion, it suffices to sliow that \± are bounded and that 
\±{x,t) are convergent as t tends to infinity for a. a. a; € fi. We recah that ip is nonnegative. 
Hence, in view also of the last of (|3.14p . we have that 



< A+(x,t) < A-(x,t) < \M\'Il(n) + for a.a. (x, t) e g^, (5.3) 



so that A± are bounded. On the other hand, it is clear that both A_ and are non-decreasing 
with respect to time, so that their convergence is ensured. □ 

Lemma 5.2. We have that: 

P < \/?* o,-^- Qoc and / \dt(p\ < +00 ; (5-4) 
pe L°°{0,oo;V) and ( \dtp\^ < +oo . (5.5) 

jQaa 



Proof. The first of (|5.4p is a consequence of Lemma [4.21 To prove the second one, we observe 
that the calculation in the previous proof yields: 



A+(t) < / A_(t) < c|17| for every i>0, (5.6) 
where c is the right-hand side of (|5.3p and |ri| is the measure of Q,. This clearly implies that 

X|9t/9p<+oo and / Xct^ < +oo, (5.7) 



because p is bounded from below. On the other hand, we have that 



^ K\dtp\^+(T 

dtip = -X -— , 5.8 

2VP 



whence immediately 

K\dtp\ 



Qt JQ 



\dM < i X = £ (A,(.) + A_ it) ^V^)<2c M 



for every t > (with the same meaning of c as before) ; the second of (|5.4p follows. To prove (|5.5p , 
we formally test p.l3p by dtp and integrate over (0,t). Then, we perform an integration by 
parts with respect to time, and obtain: 

\dtp\' +l^2^dt^+l^ Q I vp(i)p + /(p(0)) + 2 ^ 
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for every t > 0. By accounting for ()5.8|) . we see that the above equahty can be written as 
follows: 

K f {i~x)\dtp\'+ f (l-\Vp{t)\^ + fipit)))+2 f Vp^v^ 

JQt Jn \^ / Jn 

l\Wpo\^ + fipo)) +2 f ^/mvit)+ f Xa. (5.9) 
n \^ J Jn JQt 



Thus, all the terms on the left-hand side of ()5.9p are nonnegative. Moreover, the right-hand 
side is bounded, since p and are bounded and Xa E L^{Q oo) by the second of ()5.7|) . We 
immediately deduce that 

ii^x)\dtp\^ < +c», 



i.e., that the first of (|5.5p holds. Moreover, on recalling the first of (|5.7p . we see that the second 
of (l531) holds as well. □ 



Conclusion of the proof. Properties (|5.5p imply, in particular, that p is a bounded weakly 
continuous valued function on [0, -l-oo), due to the compact embedding V C H. Therefore, 
the set uj{p) is a non-empty compact subset of H. Actually, Lu{p) is also connected, due to the 
continuity of p from [0, +oo) to H and to a standard argument from the theory of dynamical 
systems (see, for instance, [U p. 12]). Then, the first properties of the cj-limit stated in the 
theorem follow. It remains for us to characterize uj{p). 

Let p°° e ijj{p) and let {f„} be a diverging sequence such that {p{tn)} converges to 
strongly in H. Then, we define pn and cpn by the following formulas: 

Pn{t) -.^ p{t + t„) and ip„{t) := ip{t + tn) for t > 0, 

and we consider their weak limits on a fixed bounded interval (0, T) (e.g., for T ~ 1). Our aim 
is to prove that such limits do not depend on time and furnish a steady state. First of all, we 
notice that 

n dtpnz+ V/9„-Vz+/ f'ipn)z- ipnPn^''^z^Q for cvcry Z G (5.10) 

JQt JQt JQt JQt 

as one immediately sees by integrating ()3.13p over (t„, i„ -I- T). Next, we observe that 

l|Pn|iL~(o,T;V) < ||pj|L-(o,oo;y) and |1</j,Jloo(q^) < ||(^||l^(q^) (5.11) 

for every n. Moreover, as ip{t) converges pointwise to <poo as t tends to infinity, we easily see 
that the whole sequence {<fn} converges to (poo weakly-star in L°°{0,T; H). Furthermore, dt^p 
belongs to i^(Qoo) by ()4.17p . By accounting for the second of (|5.1ip as well, we deduce that 

fn^^oo Strongly in C°([0,r];iP(r2)) for any p < +00. (5.12) 

Considering now (|5.1ip once more, wc infer that, to within subsequences, 

Pn Poo weakly in L°°(0, T; V) and dtPn strongly in L'^{0, T; H) (5.13) 

for some poo- We deduce that (see [H Sect. 8, Cor. 4]) 

Prt^Poo strongly in CO ([0,T];ii-), (5.14) 
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whence, in particular, Pn(0) — > p(0) strongly in H. Moreover, ()5.13|) also imply that p^o must 
be a constant, namely, Poo{i) = Poo(O) for every t S [0,T]. On the other hand, p„(0) p°° 
by assumption. We conclude that Poo{t) = p°° for every t € [0,T]. Furthermore, as / and 
any power are Lipschitz continuous on [p, we see that (|5.14l) implies that f{pn) said pn^^^ 

— 1/2 

converge to /(poo) and poo , respectively, in the same topology. Therefore, we can pass to the 
limit in (|5.10p and easily conclude that p^c satisfies (|3.19p - p.20p . 

Remark 5.3. We briefly show how to extend Theorem l3.4l to the case mentioned in Remark [3.5l 

Minor changes in the above proof are necessary. First of all, we just have to write f = a + f2 
in (j5.9p and the same equality holds true. As the sequel does not involve any special property 
of /, relations (|5.4p and (|5.5p are proved in the same way. Hence, by going through the 
above proof, we see that the only point to check is the following: if we define C„ by the formula 
Cn(t) ~ ({t + tn), does the uniformly bounded sequence {Cn} converge to a constant function (oc 
weakly-star in L°°{Qt), at least for a subsequence, and is the pair (P007C00) a solution to the 
stationary problem (|3.19p - (|3.20p ? To answer these questions in the positive, we notice that a 
weak-star limit Coo actually exists. Moreover, (|5.10p holds even for any z G L^{Q,T]V) (while 
so far we have written it just for test functions that are constant in time). On the other hand, 
the convergence (|5.14p ensures that C,oa G a[poo) a.e. in Qt and that we can pass to the limit 
in (|5.10p . We obtain that 

dtPooZ+ I Vpoo-Vz+/ (Coo+/2(Poo))2- / ^oo P^^^^ Z = (5.15) 

jqt Jqt Jqt Jqt 

for every z € L'^{0,T]V). In particular, — Apoo + Coo + ./2(Poo) ~ —{^oo/PooY^'^ at least in the 
sense of distribution on Qt, whence we deduce by comparison that Coo does not depend on time. 
Then, we immediately see that (|5.15p becomes p.27p . i.e., that (poo,Coo) solves p.l9p - p.20p . 



6 Appendix 

We here spend some words on the general forward Cauchy problem 

y'(i) + 2a(i)v/^=0, 2/(0) = yo, (6.1) 

for given a e L\^J[Q, +00) and yo G [0, +00). 

If 2/0 > 0, there is a unique, strictly positive, local solution, that has the form: 



^/y(t) = V^o - / ds, (6.2) 

JO 

as long as the right-hand side remains positive; needless to say, y may happen to tend to zero 
in a finite time. As to nonnegativc solutions, we note that every local solution can be extended 
to a global solution, because the nonlinearity is sublinear. A sufficient condition for uniqueness 
is that a is nonnegative, because the function [0, +cx)) 3 y t-^ a{t)^ is nondecrcasing for 
any fixed t. The unique, nonnegative, global solution is given by 

VW)= [yy'o- j^a{s)d^ forie [0,+cx3), (6.3) 
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where the symbol ( • denotes the positive part (for instance, if a = 1 and yo = I, we have 
that y{t) ~ ((1 — for every t > 0, and that y{t) vanishes for t > 1). Actually, to verify 

that (|6.3p always provides a solution is the matter of a simple computation. For a general a, 
the situation is more complicated, as we briefly explain. 

If a is negative, uniqueness is no longer guaranteed, because a (forward) Peano phenomenon 
might occur (for instance, if a = — 1 and yo = 0, the formula y{t) = {{t~ A)+)^ yields a solution 
for every A > 0). More generally, for yo > 0, if a solution y* vanishes a some point to and if 
a is negative in a right neighborhood of to, then there are infinitely many solutions beside y*. 
Therefore, whenever an a priori assumption on the sign of a is inappropriate (as is the case 
for the problem we study in our present paper), one would like to select the solution y that 
is maximal, i.e., that satisfies y(t) > z(t) for every t > and for every solution z to the same 
Cauchy problem. 

Now, it is known that in general any Cauchy problem for a differential equation of the 
form y' = g(t,y) has a unique maximal solution whenever is a rather general Caratheodory 
function (see, e.g., the extension of Peano's theorem [SI Thm. 4.1, p. 28] mentioned on p. 95 of 
the same book); the maximal solution can be constructed by taking the pointwise supremum 
of the solutions. For g as in (j6.ip . the square roots of the solutions describe the set 

§{a,yo) := [w G Wi];^[0, : w{0) = , w{t) > for every t > 0, 

w' = —a a.e. where ui > o|, (6-4) 
whence the maximal solution y is characterized by 

\/y{t) = sup{u'(t) : w e §(a,?;o)} for every t>0. (6.5) 
We note that every solution to the Cauchy problem (|6.ip is such that 

VW)^Vyo- I a*{s)ds, (6.6) 
Jo 

where 

a*{s):=a{s) if y{s) > Q and a*(s) := otherwise. (S-^) 

Generally, this is nothing more than an a posteriori reconstruction of y, because the definition of 
a* depends on y itself. But, if y is the maximal solution, an additional property of a* holds true, 
because a > a.e. where y = 0. Indeed, by maximality, y tries to become positive whenever it 
is possible. Thus, in view of (|6.7p . we have that 

a*{t)<a{t) fora.a. i>0. (6.8) 



At this point, we can come back to our problem. A comparison of (|2.3p with (|6.ip suggests 
that we take x as a parameter and set: 

with this, formulas (|6.5p - (|6.7p become (|3.6p - p.9p . Hence, (|3.6p - p.9p provide the maximal 
global solution to p.3p for a given p\ moreover, a sufRcicnt condition for uniqueness is that a 
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is nonnegative, a case when we have from ()6.3p that 

K{dtp{s)f + a{s) 







ds 



(6.9) 



If the initial datum is strictly positive, the last formula holds for small t even for a negative 
ct; in fact, as long as the right-hand side remains positive (see Remark 13. 3p , we simply have 
that 



■ ds . 



(6.10) 
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